Heterogeneities in the Glassy State 
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We study heterogeneities in a binary Lennard- Jones system below the glass transition using molec- 
ular dynamics simulations. We identify mobile and immobile particles and measure their distribution 
of vibrational amplitudes. For temperatures near the glass transition the distribution of vibrational 
amplitudes obeys scaling and compares reasonably well with a mean field theory for the amorphous 
solid state. To investigate correlations among the immobile and mobile particles we identify clus- 
ters and analyze their size and shape. For a fixed number of immobile particles we observe that 
the immobile particles cluster more strongly together as the temperature is increased which allows 
the particles to block each other more effectively, and to therefore stay immobile. For the mobile 
particles, on the other hand, the clustering is most pronounced at small temperatures, indicating 
that mobility at low temperatures can only be sustained in cooperative motion. 

PACS numbers: 02.70.Ns, 61.43.Fs, 64.70-Pf 



I. INTRODUCTION 

If a liquid is cooled and crystallization is avoided one 
obtains a supercooled liquid. Upon further cooling the 
system falls out of equilibrium and results in an amor- 
phous solid, i.e. a glass. During this transition from liq- 
uid to supercooled liquid to glass the dynamic proper- 
ties change drastically and they therefore have received 
much attention. Most studies are concerned with the 
supercooled liquid, for which the existence of dynamic 
heterogeneities has been well established [1] . We focus in 
this paper instead on the glass phase. 

The glassy state of matter is by definition strongly in- 
homogeneous. Different particles have different neigh- 
borhoods and consequently different dynamics. One tool 
to diagnose the heterogeneities are the vibrational am- 
plitudes In contrast to a perfect crystal these vary 
strongly from site to site in the glassy state. A simple 
way to demonstrate the heterogeneities is to pick out the 
5% particles which have the smallest and largest ampli- 
tudes. Their vibrational amplitudes are plotted in Fig. 1 
for a binary Lennard- Jones mixture (for details see be- 
low). The most striking observation is the spread in the 
data which strongly increases with increasing tempera- 
ture. Whereas for T = 0.15 the of the 5% fastest 
particles is approximately five times the {df) of the 5% 
slowest particles, this grows to a factor of approximately 
fifty close to the glass transition temperature. Thus, het- 
erogeneities are present in the glassy state for all temper- 
atures and become increasingly more pronounced as the 
glass transition temperature is approached from below. 

Our goal is to characterize these heterogeneities quan- 
titatively. Obviously the mean of the vibrational am- 
plitudes is not conclusive and we therefore investigate 
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instead the full distribution of vibrational amplitudes to 
characterize the local dynamics. We compute this distri- 
bution and show that, in the vicinity of the glass transi- 
tion, it approximately obeys scaling. We furthermore 
compare it to a mean field calculation for a different 
amorphous solid [17], the gel, and find surprisingly good 
agreement. 

To further characterize dynamic as well as time per- 
sistent heterogeneities we use various definitions to iden- 
tify mobile and immobile particles. We identify clusters 
of immobile and mobile particles where particles are de- 
fined to be connected if their relative distance is smaller 
than the position of the first minimum of the pair corre- 
lation function. To quantify the shape of these clusters 
we analyze their connectivity properties in analogy to co- 
valently bonded glasses. We determine the average size 
of the clusters, the average coordination number and the 
spatial extent as measured by the radius of gyration. In 
addition to the average values we also compute the distri- 
bution of coordination numbers as well as ring statistics. 

The majority of previous work on dynamic hetero- 
geneities are studies of mobile particles above the glass 
transition. They find that mobile particles form string- 
like clusters [2-6] which grow in length and size with 
decreasing temperature [7-9]. Immobile particles on the 
other hand do not form strings [3, 4]. In this paper we 
mostly study clusters of immobile particles below the glass 
transition and quantify their size and shape as described 
above. 

The paper is organized as follows. In Sec. II we give 
details of the simulation. Subsequently several different 
definitions of immobile particles are introduced and dis- 
cussed (Sec. III). In Sec. IV we present our results, first 
for the fraction of immobile particles (A), next for the vi- 
brational amplitudes (B) and finally the cluster analysis 
(C). We summarize our results and present conclusions 
in Sec. V. 
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FIG. 1: Vibrational amplitude (df) as a function of temperature 
for the 5% fastest and 5% slowest particles. See Fig. 5 for an 
enlargement of (d^)(T) for the 5% slowest A and B particles. 



II. SIMULATION DETAILS 

Our system is a binary Lennard-Jones mixture of 800 A 
and 200 B particles with the same mass. The interaction 
potential between particles i and j of type a, P £ {A,B} 
and at positions and rj is 

wo-"..,, -(fi^)'), (1) 

where r = jr^ — rj| and eaa — 1-0, cab = 1-5, ebb — 0.5, 
ctaa ~ 1-0, (Tab = 0.8 and ctbb = 0.88. We truncate and 
shift the potential at r = 2.5 CTq,^. From previous investi- 
gations [10] it is known that this system is not prone to 
crystallization and demixing. In the following we will use 
reduced units where the unit of length is (Taa, the unit of 
energy is eaa and the unit of time is -y/ ma\p^ j (48eAA) • 

We carry out molecular dynamics (MD) simulations 
using the velocity Verlet algorithm with a time step of 
0.02. The volume is kept constant at F = 9.4^ = 831 and 
we use periodic boundary conditions. Previous simula- 
tions [10] showed that computer simulations the system 
falls out of equilibrium in the vicinity of Tg = 0.435. We 
are interested in the glassy phase and hence study the 
system well below and close to this temperature in the 
range 0.15 < T < 0.43. 

As described in Sec. IX of [11] we equilibrate the sys- 
tem at constant temperature T — 3.0 to obtain ten in- 
dependent configurations which are at least 5 • 10"* time 
units apart. The temperature is kept constant by replac- 
ing the velocities of all particles by new velocities drawn 
from the corresponding Boltzmann distribution every 50 
time steps. Each of these configurations is then cooled 
linearly in time t {T{t) = Tq - -yt) from Tq = 3.0 to 
Tinit = 0.5 with 7 = 1.25 • 10"^ After an (NVT) equili- 
bration at T = 0.5 of 2 • 10^ time units, we then instan- 



taneously quench the system to the investigated temper- 
atures T ^ 0.15, 0.20, 0.25, 0.30, 0.35, 0.38, 0.40, 0.41, 0.42 
and 0.43. We anneal the system at each temperature 
with (NVT) simulations for 2000 time units and then run 
the production runs with (NVE) simulations for 10^ time 
units. Since the relaxation times r are much larger than 
the waiting time before the production runs (r « 8 • 10^ 
at T = 0.446), we expect to find aging effects as they 
have been studied in detail for the same Lennard-Jones 
system [12, 13] and as they have been found in related 
work [11, 14]. In this paper we do not aim for a sys- 
tematic investigation of aging effects. However, we have 
analyzed a second set of simulation runs with a different 
history in order to get a rough estimate of the effects of 
aging. The second set of runs start from well equilibrated 
configurations at Tjnit = 0.446. After an instantaneous 
quench to T = 0.15 the system is annealed for 2000 time 
units followed by a (NVT) run of also 2000 time units. 
The temperature is then raised to T — 0.2 again an- 
nealed with (NVT) and (NVE) runs each of 2000 time 
units and raised to T = 0.25 etc.. The production runs 
are for 10^ time units and are for each investigated tem- 
perature following the respective (NVE) simulations. For 
more details see [11, 14] 



III. DEFINITION OF IMMOBILE AND MOBILE 
PARTICLES 

We study heterogeneities below the glass transition 
and focus on how the particles are increasingly frozen in 
as the temperature is lowered. Even though the system 
is obviously very inhomogeneous (see Sec. I), it is not ob- 
vious how to quantify the characteristic features of these 
heterogeneities. One possibility is to study the extremes, 
the most mobile and the most immobile particles. In 
most of this paper we concentrate on the immobile ones; 
we study their organization into clusters and the size and 
shape of these clusters. 

Using four different definitions of immobility we are 
able to test if the resulting dynamical behavior is depen- 
dent on the precise definition of immobility. For two of 
the definitions we characterize the mobility of a particle 
by its vibrational amplitude [14] 



d?-|r.-fl|' (2) 

where the bar denotes a time average over the whole sim- 
ulation run. Our first definition of immobility (I) defines 
every particle i to be immobile for which 

rf- < rfcut (3) 

The second definition of immobility (II) [14] identifies 
as immobile particles a fixed percentage of A particles 
and separately a fixed percentage of B particles with the 
smallest dj. This allows us to similarly study the most 
mobile particles with the largest d^. 
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FIG. 2: Sketch of a particle trajectory to illustrate the definition 
of irreversible and reversible jumps. 




FIG. 3: The fraction of immobile particles Q as a function of 
temperature T using definition I for immobility. Q{T) is given 
both for A and B particles separately and for dcut = 0.2 and 
0.01. The simulations are with initial Temperature Tinjt = 



W2 
"cut 

0.5 



The third and fourth definition of immobility are based 

on recent work on single particle motion [11] . As sketched 
in Fig. 2a, we define a particle i to undergo a jump if 
its change in average position ARi is significantly larger 
than its fluctuation Ui {ARi > v^fj). We distinguish 
irreversible and reversible jumps where in the latter case 
the particle returns during the simulation run to one of its 
former average positions (for more details see [11]). Wc 
call all particles immobile which are not undergoing any 
jumps (III) or which are not irreversible jumpers (IV). 

Note that our definitions of immobile particles refer to 
the whole simulation run, because we want to investigate 
those structures which persist on the longest timescale, 
set by the length of the simulation run. 
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FIG. 4: Same as Fig. 3 but for definitions III and IV, for A and 
B particles, and for different initial temperatures. 



IV. RESULTS 

A. Fraction of Immobile Particles 

Let us first see how many particles are immobile ac- 
cording to the definitions I, III and IV. Fig. 3 and 4 show 
that the fraction of immobile particles 

Q = number of immobile part. / number of part. (4) 

approaches one for very low temperatures and drops 
significantly for temperatures near the glass transition. 
Q{T) depends quantitatively on the definition of immo- 
bile: the larger d^^^ the steeper is the drop of Q. The 
smallest cutoff d^^^ = 0.01 corresponds to the Lindemann 
criterion and gives rise to significant deviations in the re- 
sults for Q{T) as compared to the other definitions of 
immobility. Apparently the frozen structure can be sus- 
tained with a considerable fraction of particles whose vi- 
brational amplitude is larger than suggested by the Lin- 
demann criterion, e.g. roughly 25% of the A-particles at 
T — 0.3. This is consistent with the broad scatter of 
vibrational amplitudes in the glassy phase that we shall 
discuss below. 

The qualitative behavior is the same for all definitions 
of immobility and in particular the temperature where Q 
goes to zero is rather insensitive to the details of the def- 
inition of immobile particles (unless the cutoff is chosen 
too small, see above). Hence Q serves as a good indica- 
tor for Tg. Similarly the history of tlic^ simulation runs 
has a quantitative influence on Q{T): larger initial tem- 
perature (Tinit = 0.5) results in fewer immobile particles, 
consistent with Fig. 22 of [11]. The qualitative behavior, 
however, is not only independent of the precise definition 
of immobility but also independent of the history of the 
simulations runs. 
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FIG. 5: Same as Fig. 1 for immobile particles using definitions I 
- IV and for A and B particles. The inset shows an enlargement 
for low temperatures with linear fits through the origin. 



B. Vibrational Amplitude 

One possibility to extract the various length scales 
of the glassy state are the vibrational amplitudes {df}. 
Fig. 1 shows (c?f ) as a function of temperature T for the 
5% particles with largest (mobile, filled symbols) and 
the 5% particles with smallest (immobile, open sym- 
bols). Here (•) denotes an average over the subset of 
particles under consideration (mobile or immobile) and 
over initial configurations [15]. The remarkable spread of 
(d?) in Fig. 1 illustrates the increasing importance of het- 
erogeneities when the glass transition temperature is ap- 
proached from below (see Sec. I). In the inset of Fig. 5 we 
have included linear fits through the origin, which would 
correspond to the dynamics of a harmonic solid. We find 
for the 5% immobile particles {df){T) for T < 0.25 good 
agreement with this line fit. For larger temperatures, 
however, the immobile particles display anharmonic dy- 
namics. On the other hand, the immobile particles of 
definitions I, III, and IV include at low temperatures al- 
most all particles (see Fig. 3 and Fig. 4), and the linear 
fits through the origin are therefore only good approxi- 
mations for T < 0.2. 

The (df) of immobile particles using definitions I, III 
and IV are slightly larger but comparable to the (df) us- 
ing definition II (see Fig. 5). We find for all definitions of 
immobility the same temperature dependence. Similarly 
the analogous figures for Ti„it = 0.446 are qualitatively 
the same. 

To investigate further the heterogeneities we show 
in Fig. 6 the distribution of the vibrational amplitude 
normalized by its average. The distribution of vibra- 
tional amplitudes has been computed analytically for the 
amorphous solid which results from crosslinking macro- 
molecules, i.e. the gel state [16, 17]. In that context Q 
corresponds to the gel fraction and (df) to the localiza- 
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FIG. 6: Distribution of the normalized vibrational amplitude 
compared to a mean field theory [16]. For the simulation data 
the definition III for immobility has been used and the distribution 
is shown for A particles. 

tion length squared. It was shown that the distribution 
of localization length obeys scaling in the critical regime 
of the transition from the fluid to the amorphous solid 
state which is a second order equilibrium phase transi- 
tion. Furthermore, the scaling function has been com- 
puted within a mean field approach. The analytical re- 
sults are in good agreement with results from numerical 
simulations of crosslinked macromolecular melts [18]. 

The two systems, binary Lennard- Jones glasses and 
crosslinked polymeric melts, have completely different 
microscopic dynamics and the nature of the glass tran- 
sition, if it exists, is not understood in contrast to the 
gelation transition. Nevertheless it is an interesting open 
question to see whether the distribution of localization 
lengths P{d^) obeys scaling close to the glass transi- 
tion and how it compares to the universal function of 
the gelation transition. Fig. 6 illustrates that indeed 
for temperatures T > 0.38 P{d^) depends on tempera- 
ture only through the average squared localization length 
and that the normalized distribution compares reason- 
ably well with the universal function of the gelation tran- 
sition, as predicted by mean field theory [16, 17]. This 
agreement with theory is independent of the definition 
of immobility and independent of the history of initial 
configurations. 

C. Cluster Analysis 

The vibrational amplitude has given us an estimate for 
the length scale of single particle motion. In the rest of 
the paper we investigate how the immobile and mobile 
particles are spatially correlated with each other and how 
they are organized in clusters. 

To do so we use definitions I - IV of immobility (see 
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FIG. 7: Average number of cluster members s as a function of 
temperature. Shown are an average of s over only the biggest 
cluster of each independent configuration (sbci) and an average 
over all clusters (s) are shown for definitions I, III and IV. 

Sec. Ill) to select a subset of particles of the system. 
We stress again that immobility defined in this way is a 
time persistent feature (on the timescale of the simula- 
tion run). We then choose a particular instant of time — 
the beginning of each production run — to identify near- 
est neighbor connections among this subset of particles: 
particle j is defined to be a neighbor of particle i if their 
distance |rjj| = |rj — r^l is smaller than the position of 
the first minimum rmin of the corresponding radial pair 
distribution function of the complete system {r^in = 1.4 
for AA, 1.2 for AB and 1.07 for BB independent of tem- 
perature) . 

This definition of nearest neighbor connections gives 
rise to K distinct clusters of immobile particles, with 
two particles belonging to the same cluster if they arc 
connected. The clusters are numbered by A; = 1,2,. ..K 
and we denote by A4 the k-th cluster with Nk particles. 
The clusters of mobile particles arc defined in complete 
analogy. However, in this case the identity of the clus- 
ters refers to a particular snapshot — usually the initial 
instant of time. In contrast the clusters of immobile par- 
ticles are likely to be time persistent structures. 

To get a first idea of the organization of immobile par- 
ticles into clusters we simply determine 

{^) = {^jlN,) (5) 
fe=i 

i.e. the total number of immobile particles, Nimmob = 
J2k=i divided by the total number of clusters K, av- 
eraged over initial conditions. Here and in the following 
(•) denotes an average over initial conditions. In Fig. 7 we 
show (s) (open symbols) for the definitions I, III & IV to- 
gether with (sbci) (filled symbols), where Sbci = A^A;,bci is 
the number of particles in the largest cluster. At the low- 
est temperatures almost all particles are immobile (see 
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FIG. 8: Same as Fig. 7 but for definition II and only for (sbcl)- 



Figs. 3 and 4) and connected, so that there is only one 
large cluster and (s) = (sbci) is approximately given by 
the total number of particles. As the temperature in- 
creases the fraction of immobile particles decreases and 
the number of clusters decreases simultaneously so that 
(sbci) > (s)- Definition IV ("no irrev. jump") gives rise 
to a larger (s) for immobile particles than definition III 
("no jump"), since IV excludes fewer particles from the 
subset than III, namely only irreversible instead of also 
reversible jumpers. 

To disentangle the effects of spatial rearrangements 
and changing number of immobile particles, we use defi- 
nition II with a fixed number of immobile particles. Fig. 8 
shows (sbci) for the 5% most mobile and 5% most immo- 
bile particles as a function of temperature. In the tem- 
perature range 0.25 <T < 0.4 the 5% immobile particles 
cluster increasingly more together with increasing tem- 
perature. We interpret this such that immobile particles 
manage to stay immobile at larger temperatures by form- 
ing each others cages most successfully by clustering to- 
gether. The 5% most mobile particles on the other hand 
cluster more together at lower temperatures. A more co- 
operative motion allows the particles to be mobile at low 
temperatures. 

This picture gets further support from the observa- 
tion that with increasing temperature the number of clus- 
ters {K) increases for the 5% most mobile particles and 
decreases for the 5% most immobile ones (see Fig. 9). 
This temperature dependence is qualitatively the same 
for simulations with different histories (Tjnit = 0.446 and 
Ti„it = 0.5). 

Another interesting question concerns the percolation 
properties of clusters of either immobile or mobile par- 
ticles. We vary the percentage p of immobile particles 
and thereby the clusters Afk as well as the total number 
of clusters K. For a fixed temperature T we compute 
the fraction Pporc of the ten independent configurations 
with a percolating cluster of immobile particles. We de- 



6 



25- 



20- 



^"15- 



10- 



^ , I , , ^ 

O ^ r. . =0.446 5% immob. 

init 

A- -A T.^^=0.446 5% mobile 

»r. . =0.5 5% immob. 
init 

1 T. =0.5 5% mobile 




0.2 



0.3 



0.4 



FIG. 9: Average number of clusters {K} as a function of tem- 
perature using definition II. 
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FIG. 10: The average cluster size (5') as a function of percentage 
p using definition II for immobility. The average is over all non- 
percolating clusters and over independent initial configurations. 
The inset shows the fraction of percolating configurations as a 
function of p. 



fine a cluster to be percolating if its maximum extension 
= max \ri^ — ^jj between any pair of particles i and 
j satisfies ri, > 8.0 for any direction u G {x,y,z} [19]. 
As can be seen in the inset of Fig. 10 this fraction rises 
sharply from zero at a critical value Pc, which is approx- 
imately independent of temperature. A central quantity 
from percolation theory is the average cluster size, de- 
fined as 



{S){P) 



K 



with ns = '^5Nk,s (6) 



Percolation theory predicts a divergence of {S{p)), as the 
percolation transition is approached from either side (ex- 



'a 



-1- 



J , , , I L. 



T .=0.5 

imt 




o 


T=0.15 


□ 


T=0.20 


o 


T=0.25 


A 


T=0.30 


V 


T=0.35 


> 


T=0.38 




line fit 



FIG. 11: Log-log plot of the radius of gyration Rq as a func- 
tion of the number of cluster members s using definition I for 
immobility. The line is a linear fit with slope 0.55. 



eluding the percolating cluster for p > pc). In our case 
we observe a strong increase of (S), a maximum around 
Pc and a subsequent decrease for p > Pc (see Fig. 10). 
The strong increase of {S) around pc is reminiscent of a 
percolation transition with, however, a threshold value 
that is significantly lower than the corresponding value 
for hard spheres [20] . This may be due to the strong cor- 
relations between immobile particles in addition to the 
correlations due to the interaction potential V{r). Fur- 
thermore our analysis is severely hampered by the small 
size of our sample, which percolates much easier than a 
larger sample, so that pc is underestimated. 

The spatial extent of clusters of immobile and mobile 
particles can be quantified with help of the radius of gy- 
ration. The latter is defined for a particular cluster Mk'- 



with the center of mass of the cluster 

1 



(7) 



(8) 



To study the scaling of Rq as a function of the number 
of particles in a cluster s, we show in Fig. 11 a log- log 
plot of RoiAfk) versus s{Afk) for all clusters A/fe of all 
independent configurations. We find that approximately 
Rq oc s^/'^-f with a fractal dimension df « 1.8 (and sim- 
ilarly for the definitions II - IV). If the clusters were 
approximately straight lines one would expect df = 1.0, 
if they were compact df = 3.0. Our result points to frac- 
tal structures, however this should be taken with caveat 
since the clusters are too small to observe true asymp- 
totic scaling behavior. 

Next, we investigate the connectivity within single 
clusters. Since two sites have been defined to be con- 
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FIG. 12: Average coordination number of members of the 
biggest cluster (^bcl) 3S a function of temperature using defi- 
nition II. 



nected if sufficiently close in space, the connectivity is 
not a purely topological feature, but gives us information 
about the shape of the clusters. A convenient measure, 
which does not require a large variety of cluster sizes, is 
the coordination number z. We define the average coor- 
dination number z{J\fk) of a cluster A/fc as 



1 



ieAfk 



(9) 



where Zi is the number of neighboring particles j € A4 of 
particle i. Fig. 12 gives a plot of (zbci) = (•^(A/'bci)), the 
average number of neighbors of particles in the biggest 
cluster averaged over independent initial configurations. 
We observe that (zbci) is increasing for the 5% immo- 
bile particles with increasing temperature, independent 
of the history of the runs. The immobile particles block 
each other at larger temperatures by clustering more to- 
gether and by forming more compact clusters. Clusters 
of mobile particles on the other hand become more elon- 
gated with increasing temperature. Although (zbci) is not 
equal to two as it would be for strings, (zbci) approaches 
2.5 as T — > Tg. This is reminiscent of the results [7-9], 
where increasingly string-like motion for mobile particles 
were observed as the glass transition is approached from 
above. (The authors of [7-9] use a slightly different defi- 
nition of mobility than in this paper and define a string 
via velocity- velocity correlations.) 

More information about the connectivity and the shape 
of the cluster can be found from the distribution of co- 
ordination numbers of the particles in the biggest clus- 
ter for each initial configuration /'(zi^bci) — ^'(•ZiGA/'bci)- 
Fig. 13 illustrates that for low temperatures z — 3 and 
z = 2 are the most common coordination numbers of the 
mobile particles in the biggest cluster. However, for high 
temperatures z = 2 dominates and the variance of the 
distribution is smaller, indicating string-like clusters. In 




FIG. 13: Distribution of coordination numbers bcl) of the 
5% fastest particles in the biggest cluster for simulations with 
Tm\t = 0.5. The inset shows P{zi [^^\) for the 5% slowest particles 
for simulations with Tipjt = 0.446. 



the case of the 5% immobile particles, on the other hand, 
the peak is at z = 4 and the variance of P{zi^j^^^^) is 
strongly increasing with increasing temperature (see in- 
set of Fig. 13). This indicates increasingly more compact 
clusters with increasing temperature. 

We finish our cluster analysis with ring statistics. 
Whereas ring analysis usually is applied to network form- 
ing glasses, it is here applied to the binary Lennard- Jones 
system, which is a fragile glass former [21]. In the fol- 
lowing we call any path a set of nearest neighbor con- 
nections. The most commonly used method, introduced 
by S. V. King [22], searches for rings by finding for each 
particle and any pair of its neighbors the shortest path 
between them. As pointed out by Franzblau [23], this 
definition has serious disadvantages such as not counting 
certain "intuitive" rings, which leads Franzblau to the 
"shortest-path" (SP) criterion: an SP ring contains for 
each pair of particles in the ring the shortest path be- 
tween the pair in the cluster under consideration. We 
use the SP criterion and search algorithm of Franzblau 
to find within the largest cluster all rings r = IjA'^.bci 
each of size rir, i.e. with Ur connections. 

Fig. 14 shows the distribution P{nr) for the biggest 
cluster of the 5% immobile and 5% mobile particles. We 
find that rings are dominantly made up of three particles, 
only a few rings are longer. This means that in the case 
of the "stringier" mobile particle clusters, the clusters are 
mostly not closed rings. 

In the case of immobile clusters using definitions I, 
III and IV larger clusters are formed leading to larger 
average ringsize 
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(10) 



(see Fig. 15). For low temperatures almost all particles 
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FIG. 14; Distribution of ring sizes P{nr) within the largest clus- 
ter of each independent configuration using definition II. 
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FIG. 15: Average ring size (ribci) as a function of temperature 
using definitions III and IV and the comparison of (n) for the 
complete system (100%). 



are members of the biggest cluster of immobile particles 
and (nbci) is therefore comparable to (n) of the whole sys- 
tem. For larger temperatures the clusters become more 
compact and on average larger ring sizes occur. 



V. CONCLUSIONS 

We investigate a binary Lennard-Jones system below 
the glass transition. To study the inhomogeneous struc- 
ture of the glassy state quantitatively wc focus on the 
most mobile and immobile particles. Since the identifi- 
cation of these is not unique — at least not in a simu- 
lation run of finite length — we test different definitions 
of immobility. The first two definitions are based on the 



vibrational amplitude and the last two define a particle 
to be immobile if it is not jumping. 

The length scales of single particle motion varies 
strongly among the particles, giving rise to a broad distri- 
bution of vibrational amplitudes for mobile and immobile 
particles. The range of vibrational amplitudes increases 
as the glass transition temperature is approached from 
below: The system does not only become more heteroge- 
neous with increasing temperature, but also the range of 
lengthscales extends to much larger scales in the vicinity 
of the glass transition than at low temperatures. Fo- 
cussing on the immobile particles only, we observe the 
distribution of vibrational amplitudes to follow a scaling 
plot close to the glass transition. It even compares rea- 
sonably well with a mean field theory for the amorphous 
solid state, developed in the context of the gelation tran- 
sition. Tentatively this may be taken as evidence for the 
recently suggested picture of the more general jamming 
transition [24] , including both the gel transition and the 
glass transition. 

To characterize the spatial arrangement and correla- 
tions among immobile and mobile particles, we identify 
clusters of particles, which are defined to be connected if 
their relative distance is smaller than a threshold value, 
given by the first minimum of the pair correlation func- 
tion. Whereas for the immobile particles these clusters 
are approximately time persistent and hence characterize 
the static heterogeneity, the clusters of mobile particles 
are relevant for the dynamic heterogeneity of the glassy 
state. 



If immobility is defined by the vibrational amplitudes 
or the jump characteristics, then the number of immobile 
particles increases strongly as the temperature is lowered. 
This effect dominates the statistics of cluster number, size 
and shape. Baljon et al. [25] find similar results for thin 
polymer films by looking at the percolation probability 
of immobile monomers as a function of temperature. It is 
instructive to keep the number of immobile/mobile par- 
ticles constant, namely the p% particles with the small- 
est/largest vibrational amplitude. The p% mobile parti- 
cles on the other hand cluster more together at low tem- 
peratures and form more elongated string-like clusters 
at larger temperatures. Similar strings have been found 
for the same system above the glass transition [2, 3] and 
for a similar system above and slightly below the critical 
temperature [5]. These simulations, our results as well 
as experiments [9] point to string-like clusters of mobile 
particles which grow close to the glass transition tem- 
perature, both when approached from above as well as 
when approached from below the glass transition. We 
find that with increasing temperature the p% immobile 
particles cluster more strongly together and form more 
compact clusters which allows the particles to block each 
other more effectively, and to therefore stay immobile. 
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